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Using the Bak-Sneppen model of biological evolution as our paradigm, we investigate in which cases noise can be 
substituted with a deterministic signal without destroying Self-Organized Criticality (SOC). If the deterministic 
signal is chaotic the universality class is preserved; some non-universal features, such as the threshold, depend 
on the time correlation of the signal. We also show that, if the signal introduced is periodic, SOC is preserved 
but in a different universality class, as long as the spectrum of frequencies is broad enough. 



I. INTRODUCTION 

Due to Nature's inherent complexity, a lot of effort has gone 
ipfrO) develo ping mathematical models to describe it, if only qual- 
itatively Jl]-[l6|]- Among all natural processes, evolution has at- 
tScfcd a lot of attention because of its global consequences for 
lM^[ |r^ , |l8| ]. Within the realm of evolution, one of the most 
ferjKdly argued topics is that of the explanation of mass ex- 
tinctions |p^| . Indeed, from a gradualistic point of view, mass 
extinctions are rare events, due mainly to external abiotic fac- 
tors* such as earthquakes, meteorites etc. poj . From the point 
of view of punctuated equilibrium, on the other hand, mass ex- 
t jsqfcibns are bursts of activity between periods of stasis (2l] |2^| . 
Th^ fossil record shows that the distribution of mass extinctions 
f^Dbws a power law ]2(| . Among the many models proposed to 
describe evolution, those exhibiting Self-Organized Criticality 
(SOC) p7| , p8[ are of particular interest. In layman terms, a 
svsiem is called self-organized critical when it evolves towards 
already state in which certain physical quantities show fluc- 
tjjy^ions on any space and time scale (they follow power law 
cflg^ribution) . 

4[J particular, we will concentrate on a model for macroevolu- 
ti(£« proposed by Bak and Sneppen (BS) [p5| . In it, extinctions 
ara associated with avalanches of activity without an inherent 
time or length scale. In the original version of the model, no in- 
fliffiJice of the environment was taken into account if not implic- 
itly^n the fitnesses of the species. Later, Newman and coworkers 
introduced a modified version of the BS model in which an 
environment al stress is introduced. All these versions of the 
rti^del show self-organized criticality and the randomness in the 
microscopic rule seems to play a relevant role. 

this paper we show that, if one substitutes the random up- 
dating of the dynamic variable with a chaotic or a periodic map, 
SOC is not destroyed p9| |. Some non-universal features will, 
however, depend on the time correlation of the signal. More- 
over, if the signal introduced is periodic, SOC is preserved but 
in a different universality class, as long as the spectrum of fre- 
quencies is broad enough (a brief discussion of some of these 
results can be also found in ]3l|]). Similar conclusions have been 
drawn in the context of standard fluctuation-dissipation pro- 
cesses by Bianucci et. al. ]3^-|34|] . The paper is organized as 
follows. In Sec. |l| we review the main features of the model as 
they are presented in the literature. After introducing the maps 



in Sec. Ill, we show in Sec. IV the results obtained for different 
deterministic updating rules. Our conclusions together with a 
brief description of some open problems can be found in Sec. |v|. 

II. GENERAL DESCRIPTION OF THE MODEL 

The Bak-Sneppen model describes an ecosystem as a collec- 
tion of N species on a rf-dimensional lattice, each one of which 



can have T traits associated to it pj|. To each one of these 
traits corresponds a fitness described by a number / between 
and 1 p6| ]. Here, for simplicity, we consider the case with one 
trait, d = 1, and periodic boundary conditions. To fix nota- 
tion, we consider a one dimensional lattice of length N. The 
initial state of the system is defined by assigning to each site j 
a random fitness /g chosen from a uniform distribution. The 
dynamics proceeds in three basic steps: 

1. Find the site with the absolute minimum fitness on the 
lattice (this site will be called the active site) and its two 
nearest neighbors. 

2. Change, at the same time, the values of their fitnesses by 
assigning to them new random numbers from a uniform 
distribution. 

3. Go to step 1. 

After an initial transient that will be of no interest to us, a non- 
trivial critical state is reached. This critical state, characterized 
by its statistical properties, can be understood as the fluctuat- 
ing balance between two competing "forces" . Indeed, while the 
random assignation of the values, together with the coupling, 
acts as an entropic disorder, the choice of the minimum acts 
as an ordering force. As a result of this competition, at the 
stationary state the majority of the f J have values above a cer- 
tain threshold f c . Only a few will be below f c , namely those 
belonging to the running avalanche (see |^5|,^9) for a detailed 
discussion). Since the avalanches are the basic and fundamental 
mechanism of the model it is therefore worth describing them 
in more general terms. Let us suppose that the system is al- 
ready at stationarity and let us find the minimum fitness, say 
jo < fc- We update it together with its two nearest neighbors 
(the actual value of the minimum does not really matter). This 
updating creates disorder in a small region in space, where most 
probably there are some lattice sites with P < f c . Then, the 
new minimum will most probably be among the last three sites 
changed. The active site most probably will be one of the two 
nearest neighbors, thus affecting another site to the left or to the 
right. In the following time-steps new sites will be touched by 
the avalanche of mutations. Here one sees clearly the two afore- 
mentioned forces at work: Disorder (every new value is chosen 
at random) and order (we decide to mutate always the smallest). 
Since the equilibrium drives the threshold to the value f c > f a , 
this means that the / -avalanche (an avalanche during which 
all the selected minima are below f a ) will eventually come to 
an end, in a finite time s/ o = s Q . During this process the f a - 
avalanche will also cover a certain number of lattice sites, i.e. 
it will also have a spatial size n[° v = n . This feature gives the 
possibility to analyze and find the critical values even without 
considering the system or lattice as a whole but simply analyz- 
ing the statistics of / Q -avalanches fjofl . Moreover, the avalanche 
dynamics shows that, as long as f Q is close to f c the average size 



1 



(both in time and space) of the related avalanches has to grow 
considerably. Both averages will eventually become infinite as 
fo — ¥ fci but this does not mean that all the avalanches are 
infinite (or of the maximal space or time lengths allowed by the 
simulations). These / -avalanches can be described by means 
of a distribution function |E3,B9| 



Pa(sJo) = S- T F{sUc-fo) 1/a )- 



(1) 



In Eq. (0) s is the time-size of the avalanche and F is some yet 
unknown scaling function that behaves like 



F{sUc~fo) 1/a ) = 



1 as s — 

as s>_p = (f c 







The average duration of an / G -avalanche is given by 

So = (fc ~ fo)- 1 , 



(2) 



(3) 



where the exponent 7 is given in terms of the previously defined 
exponents r and a by 



7 = 



(4) 



Numerical calculations provide good estimates for the value of 
the threshold as well as the two exponents r and a |2^,^] (see 
Table 1). It is also useful to define other exponents that can be 
easily obtained from numerical simulations. First we consider 
the first return time distribution Pf(n), namely the distribution 
of the times between two consecutive updatings of the same site 
(when it is the minimum). Another distribution function is the 
all return time distribution P a u (n) , namely the probability that 
a given site, active at time t = 0, is active again at time t. In 
both cases, one defines the corresponding exponents by 



P f (n) ~ rT T ' and P aU (n) 



(5) 



In Table 1 we have listed the values of these exponents as 
they are given in the literature. It is worth noticing that these 
exponents are not independent quantities. Indeed, the scaling 
relations derived in p5|,p9[ show that at most two of them can 
be independent. However, using the master equation 



(Wo) 



dPa(sJo) 
dfo 



-Pa(sJ )n (s) (6) 



+ ^2 -^(si, fo)n (si)P a (s - Si, f ) 

S1 =l 

derived in Efl and the fact that n (s) ~ s Ta " and a = 1 + 
T a — t p5|,39|, one proves that the only independent exponent 
of the model is T a u of Eq. (||) Ji2| . From Eq. (Q) one can also 
derive an infinite hierarchy of equations for the moments of the 
distribution. The first equation in this hierarchy 



d log s 
dfo 



l-/o 



(7) 



relates the exponents (||) to the average number of covered sites; 
here n is the non-universal average number of sites covered by 
the / -avalanches. Putting Eq. (||) into Eq. (Fij) gives the so 
called 7-equation pa,p9[, 



7 = 



l-/o 



(fc - fo) ■ 



(8) 



For models belonging to the same universality class i.e. with 
the same 7 this equation relates the non-universal quantities n 
and the threshold f c . In particular, as we shall see in Sec. IV A, 
to a bigger f c must correspond a smaller h for fixed f a . 

An interesting consequence of Eq. (|sj> is that it is possible 
to change f c while remaining in the same universality class. 
This can be obtained by modifying the entropic tendency. In- 
deed substituting the random updating with a correlated chaotic 
system one introduces a correlation that leads to an increase 
towards 1 of the threshold. On the other hand, a greater corre- 
lation in the updating map means that the system spends more 
time in the same site, thus covering less sites in the same number 
of time-steps in comparison with a less correlated map. From 
Eq. (pi) it is clear that as n decreases, f c increases. 



quantity 


value 


error 


fc 


0.66702 


8 


T 


1.073 


3 


a 


0.343 


4 


7 


2.70 


2 


T f 


1.58 


2 


Tall 


0.42 


2 



Table 1 The first four values exponents are quoted from 
while the last two from (25U33]. 



III. MAPS 

As we have seen before, the source of mutations in the Bak- 
Sneppen model is the presence of random noise in the system. 
Since a chaotic map may exhibit statistical properties similar 
to those of random noise, a similar competition between order 
and disorder could be established when one substitutes random 
updating with chaotic updating. To understand the similarities 
as well as the differences between the two kinds of updating, in 
this section we discuss some general properties of maps. 

A deterministic map is a rule in which the new value of the 
variable is given by 



fi+i = nil) 



(9) 



with F a deterministic function and j the lattice site. In what 
follows, we will only consider maps of the unit interval onto 
itself (usually called unimodular maps). Disregarding periodic 
trajectories, one can define several statistical quantities that are 
generally used to describe the properties of a generic sequence 
{/«}• 

The first quantity of interest to us is the invariant measure, 
fj.(f). Formally, the invariant measure for a unimodular map is 
defined by 



1 N 

H f0 (f)= Hm **(/„)] 

2=0 



(10) 



If Hf (f) does not depend on the initial value fo, the map is 
called ergodic (and one refers to the measure as fi(f)). If a 
system is ergodic, time averages are equivalent to phase space 
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averages, and then the time average of any function g(f) can be 
computed as a phase space average via 



< g 



>= lim 

N — >cx 



N 
i=0 



g(fi) = / rifMf)df. 



(ii) 



To describe the behavior of individual trajectories one needs 
more detailed information provided by the Lyapunov exponent 
A. The Lyapunov exponent measures the average rate of sepa- 
ration in /-space of two given trajectories per unit of time. It 
can be computed as 



A = J im ^Ekel^WI 

jv^oo W at 

i=0 J 



(12) 



If a map has a Lyapunov exponent A > 0, this means that two 
trajectories will diverge from each other exponentially. In this 
case the map is called chaotic. This property has a very impor- 
tant consequence: A very small perturbation in the initial con- 
dition will produce a completely different outcome. Moreover, 
successive outcomes of a chaotic map will behave like a stochas- 
tic variable (statistically speaking). Finally, we will make use 
of the autocorrelation function C (m) , defined as 



C(m) = ^(/, 



/)(/» - /) ■ 



where 



/ = 



1 N 



(13) 



(14) 



»=o 



The function C (m) is a measure of how the deviations from the 
average at time i are related to the deviations from the average 
m steps apart [ff~I|| . In particular, chaotic maps are expected to 
show exponentially decaying autocorrelation functions, i.e. 



C(m) 



- TI L J T 



(15) 



where r is the correlation time. 

We will now proceed to summarise the properties of the dif- 
ferent maps we will be using in the Bak-Sneppen model. Before 
continuing, it is worth mentioning that in principle, the case 
of random noise can be considered as a particular case of 
in which F(f n ) = ip(n) with ip(n) a random variable with a 
uniform probability density J43| . 



A. Bernoulli maps 

Let us start by considering the Bernoulli map 
fi + i = G r (fi) = [rfi], 



namely 
(16) 



where [/] stands for the value of / modulus 1 and r G N, r > 1 is 
a constant. It has been shown (see H] and References therein) 
that this map has a uniform invariant measure 



(BM) 



(/) 



1 . 



(17) 



where the function fi r (f) has been defined in Eq. ([To|). Moreover 
this map is chaotic and is characterized by a Lyapunov exponent 
given by 



A< SM ) =logr. 



(18) 



For this map one can easily compute the time autocorrelation 
function, namely 



1 



12r r ' 



1 

12 f 



-m In r 



(19) 



where C tyBM \m) has been defined in Eq. 
tion time is given by 

1 

lnr 



3f) and the correla- 



(20) 



One sees that the correlation time decreases as r increases. This 
means that given two maps with different values of r, the one 
with the bigger value of r will be closer to true random noise 
and then will de-correlate faster. As we shall see in the following 
section, this last property is of crucial importance in order to 
understand the differences between BS models with different 
Bernoulli maps. 



B. Logistic map 

Let us now consider the logistic map (sometimes called 
Feigenbaum map), namely 



/n+l 



(21) 



The reasons for studying this map are manifold. On the one 
hand, this map has already been considered in the context of 
biological evolution models and population dynamics [|^]-[l7],|o| 
and can thus provide a possible deterministic interpretation of 
the evolution inside every ecological niche. Moreover, it has 
been shown that it describes the behavior of a wide variety of 
systems in nature E3]. On the other hand, it has a regime in 
which it is chaotic as well as one in which it is not, depending on 
whether A is bigger or less than the critical value Aoo ~ 3.56994 
|44| (for A > Aqo there are windows in which the map is periodic; 
in this paper we will take A outside these windows). 

If we consider the particular case A — 4, the invariant density 
for this map is given by 



n {LM Xf) 



1 



- f) 

the Lyapunov exponent is 

A (LM) =ln2 

and the correlation function is given by 



(22) 



(23) 



(24) 



The fact that this map is chaotic does not mean that the tra- 
jectory cannot be written explicitly. Indeed, it is easy to see 
that 



f n = sin 2 (2"c) 



(25) 



with the initial condition /o = sin 2 (c), is a trajectory of the 
logistic map in the case of A = 4. 
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C. Tent Map 



sin(w,n + c6,-) + 1 



(30) 



To better illustrate the effects of time correlations in the up- 
dating, we will also need the so-called "tent" map (a "linear 
version" of the logistic map), defined as 



Jn+l 



p < - 



2(1 - P n ) Si > | 



(26) 



This map is chaotic with Lyapunov exponent A™ = In 2. Con- 
trary to the case of the logistic map, the invariant measure for 
the Tent map is uniform i.e. 



M (TM) (/) = 1, 
with an autocorrelation function given by 

like in the case of purely random noise. 
1.0 



F(f) 



(27) 
(28) 




0.0 



FIG. 1. Shifted Tent map corresponding ton = 1/3. In particular, 
for 77 = one obtains usual Tent map, Eq. (E6I). 



For our applications (see next section), we find useful to define 
a modified version of the Tent map (Eq) ("shifted Tent map") 
in which we cut the y — axis at 1 — 77 and then shift the function 
rj up modulus one, as shown in the following equation 



fn+1 



1 km < p < I 



77+1-2/^ l<fl<ktn 
[ r, + 2-2fl l±2</i<i. 



(29) 



As an example, in Fig. [I] we can see the plot of Eq. ( p9| ) for the 



case 77 



D. Periodic Map 

So far we have only considered chaotic systems. However, 
there are maps that are not chaotic but are ergodic. Let us 
consider a simple example of such a case in which the "signal" 
is provided by an integrable system, that is a sequence given by 



where the <jj, (j>j's are the angular frequencies and initial phases 
respectively. This can be rewritten as a map of j 3 n onto f^+i as 



j _ sin(arcsin(2/^ — 1) + Wj) + 1 

Jn+l ~ 2 



(31) 



where the initial condition /q is given by 



j _ 1 + sin (j)j 



ft 



The invariant measure is not constant, it is symmetric around 
/ = 1/2 and peaked close to / = 0, 1 namely 



Mf) = 



1 



(32) 



Since this "signal" is not chaotic, the Lyapunov exponent is zero 
and the correlations will not decay exponentially. In fact, the 
correlations are given by 



(33) 



At this point it is worth emphasising that these correlations are 
correlations for a given sequence. If we consider two sequences 
with different values of u), <f> the correlation will be different. 



IV. MODELS 

Through the use of the different maps presented in the previ- 
ous section, we shall show here that the random updating is no 
longer a necessary requirement to have SOC. Moreover we will 
also show that as long as the map at hand is chaotic the system 
does not change the universality class, i.e. all the exponents are 
the same as in the case of random updating. 

While the presence of critical behavior was somehow ex- 
pected, it is still surprising that the universality class does not 
change. This means that the system is able to self-organize at a 
higher level: It takes into account the temporal correlation (or 
the average time spent in every site) by increasing the threshold, 
so as to have the same statistical properties. What is even more 
remarkable then is the fact that close enough to the threshold 
it is not possible to distinguish the random updating case from 
the chaotic one from the microscopical point of view, being the 
statistical properties and all the variables exactly the same. As 
a consequence, all the equations and relations shown in Sec. || 
are still valid for all the cases with chaotic updating. 

We will show this equivalence through an infinite sequence 
of models with Bernoulli updating, logistic and Tent map up- 
dating. In fact, the same kind of analysis performed on the 
case with the (modified) Tent map can show that the time cor- 
relations are actually the ones responsible for the shifts in the 
thresholds. 

However, the universality class is not always preserved. In 
fact, if one chooses a non-chaotic updating rule the critical ex- 
ponents may change. We will show that by considering quasi- 
periodic updating rules (pl|). 
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A. Bernoulli updating. 

Let us consider a chaotic updating rule, whose statistical 
properties resemble those of a stochastic function, namely the 
Bernoulli map, Eq. (fl6|). 

In Fig. ^ we show the power-law behavior of the first and all 
return probability distributions in the case r = 2. The critical 
exponents obtained coincide with those found in |^5|,[39| for the 
random updating. Moreover, for all values of r the system falls 
in the BS universality class, i.e. it always has the same critical 
exponents. 



SOC is still preserved within the BS universality class. How- 
ever, in this case, the distribution of the generated numbers is 
not uniform and consequently it influences the distribution of 
the fitnesses at the stationary state. 

Turning now to Fig. ||, we can see that the thresholds for 
the Bernoulli updating approach the BS value as r increases. In 
Fig. |I] we have also plotted the best fit we could find for the 
curve (f c (r) — fc S )- This fit, that corresponds to a power law 
r~ ' 78 , still remains an open problem from the theoretical point 
of view. 
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Time t 



FIG. 2. First and all return distributions (non normalized) for 
a BS model with Bernoulli updating rule with r = 2. For all the 
simulations shown here, we used a lattice of 2 14 sites and 5 x 10 9 
iterations exploiting the tree-algorithm explained in [35]. 
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FIG. 3. Distribution of the fitnesses for r = 2, 3, 7, 10; the thresh- 
old for r — 2 is quite different from the usual BS threshold while 
the threshold corresponding to r = 10 is very close to the BS value 
(given by the vertical line). For all the simulations shown here, we 
used a lattice of 2 14 sites and 5 x 10 9 iterations. 

The stationary distribution of the fitnesses, on the other hand, 
follows a different pattern. Indeed, Fig. || shows that the thresh- 
old for r = 2 is bigger than the one found for the random case. 
On increasing the value of r, the threshold moves towards the 
BS value (see Fig. ph. For non integer values of r (r > 1), 
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FIG. 4. Values of the thresholds f r for r — 2, . . . , 31 after sub- 
tracting the value 0.656 for the random updating (Th — f r — /as). 
For all the simulations shown here, we used a lattice of 1000 sites 
and 5 x 10 4 measurements. 

There is a qualitative explanation for this behavior of the 
thresholds. As we briefly mentioned in Sec. 0, the change in 
threshold is an indication of the correlation in time of the map 
we are using. Indeed, by looking at Eq. (||) we can see that an 
increase in f c corresponds to a lower value of n Q for fixed f . 
This fact means that the system spends more time per site and 
this reflects the fact that it needs more time to de-correlate. At 
this stage one can also ask if Eq. (^J) remains valid even with 
a correlated map and if it is not necessary to introduce non- 
universal factors. The answer is given by noticing that for fixed 
distance A/ from the threshold the value fiA f is the same for 
all models, leading to 



n Af 



1 

A/ 



(34) 



where 7 is the universal exponent introduced in Eq. (0) | pl[ . 
This means that looking at the system from a distance A / from 
the threshold it is not possible to distinguish two systems which 
have the same critical exponents. 



B. Logistic and Tent updating. 

In previous subsection we showed that models with time- 
correlated updating self-organize into a stable configuration 
with a threshold bigger than the one in the random updating. 
At this stage, it is natural to consider updating that, even if 
deterministic, are <5-correlated. In particular, we consider up- 
dating rules given by the logistic and Tent maps. 
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FIG. 5. First and all return distributions (non normalized) for a 
BS model with the logistic map with A = 4 as updating rule. The 
exponents are the same as for the BS model Tf — 1.58 and r a = 0.42. 
In all the simulations shown in this figure, we used a lattice of 2 13 
sites and 10 s iterations. 
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FIG. 6. Distribution of the fitnesses for the logistic map. One can 
easily see the effect of the non-uniform invariant measure. 

Let us start by taking as updating rule for the fitnesses the 
logistic map, Eq. As Fig. [| shows, for those values of A 

for which the map is chaotic, the system not only exhibits SOC 
but also stays in the same universality class as the original BS 
model. For A < Aoo we find that the system is not critical any 
more. This is due to the fact that, for A < Aoo the map goes to 
a periodic orbit, and consequently the updating is not ergodic. 
This case is equivalent to a BS model with finite number of 
states for the fitnesses. In terms of our previous picture, the 
disorder force is too weak to ensure SOC. 



FIG. 7. The main body of the figure shows the distribution of the 
fitnesses for the Tent map updating Eq. (|2^). The plot in the inset 
shows the fitness distribution for the case of the shifted Tent map. 

One characteristic of the logistic updating is that, since the 
invariant measure is not uniform, the distribution of the fit- 
nesses above threshold is not uniform (see Fig. ^). This is not 
the case for the Tent map, Eq. (pq ) or the shifted version of 
it, Eq. (p9|). For both these cases the fitness distribution in 
the critical state is shown in Fig. fj]. One observes that there 
is a peak in the fitness distribution in the neighborhood of the 
threshold for Eq. (p6[). This can be understood as produced by 
the interplay between the dynamics of the updating rule and the 
Bak-Sneppen dynamics. Indeed, the Tent map has an unstable 
fixed point at /* = | > jf s . Then if a site right below thresh- 
old is chosen as the minimum, the updated value will be above 
threshold, but still close to it. The next update in the same 
site will put the value of the fitness again below threshold, if a 
little bit further apart Q). Then, one needs to update this site 
several times to remove it from the neighborhood of the thresh- 
old. Consequently, the probability of finding a site with fitness 
in the neighborhood of f c is higher than in the random update. 
If, on the other hand, one introduces the shifted Tent map, Eq. 
(p9|), where the fixed point is not close to f^ s , the distribution 
of the fitnesses above threshold is uniform, resembling exactly 
the random updating case (see Fig. 0) . 

Comparing the different chaotic maps we can draw several 
conclusions. First, time correlations in the updating immedi- 
ately reflects in a shift of the threshold in the sense that to 
higher correlations correspond higher thresholds. Second, as 
shown by the shifted Tent map, the other higher correlations do 
not in principle produce any measurable change in the statistical 
properties of the system. 



C. Periodic updating. 

Since time correlations in the updating rule do not, in princi- 
ple, destroy SOC, it is worth considering systems in which the 
time correlation of the updating does not decay exponentially. 
As shown in pi}] , the simplest example of this class is given by 
a model in which the choice of the new fitness is done according 
to the periodic map, E q. (|3l| ). 

As mentioned in Sec. |II]|, choosing the initial phases is equiv- 
alent to choosing the initial condition of the system. Conse- 
quently, we take our phases <f»j at random (0 < <fij < J). Our 
simulations indicate that if the frequencies are the same, that 
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is uij = u>, the strong synchronization of the sites along the 
lattice destroys criticality (even though the fitnesses are orga- 
nized above a threshold). Indeed, the system develops a typical 
scale that is observed in the way of a cutoff in the distribution 
probabilities. If we now choose the frequencies ujj such that 

ujj ^ uj 1 (35) 

the situation changes dramatically. If we characterize the fre- 
quency distribution by two numbers, namely its centre luq and 
its width Aoj, the behavior will indeed depend on both. 
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FIG. 8. (a) First and all return distributions (non normalized) for 
a BS model with disordered periodic updating rules; the exponents 
are r/ = 1.65(1) and r a = 0.38(1). (b) Distribution of the fitnesses. 
In all the simulations shown in this picture, we used a lattice of 2 12 
sites and 5 x 10 9 iterations. 

Even after long numerical investigation, the exact functional 
form of this dependence cannot be outlined in a satisfactory 
way. Nevertheless, it is clear that over a whole range of values of 
the two parameters the system recovers a critical behavior. An 
example can be seen in Fig. [§, where we chose luq = 19.57T and 
Au) = 19.57T. As mentioned above, the universality class changes 
with respect to the original BS model, with tj — 1.65(1) and 
r a = 0.38(1), but the SOC behavior is preserved. We observed 
that this universality class depends on the values of u>o and Auj. 
For the sake of clarity, we show here only one example, out of 
many, that illustrates the point. 

V. CONCLUSIONS 

Self-organized criticality in the BS model comes from the com- 
petition between the disorder in the updating and the ordering 
pressure of the minimum rule. For a given lattice and a given 
set of dynamical rules, the use of stochastic updating is tanta- 
mount to the introduction of maximum disorder. On the other 
hand, chaotic maps produce series of numbers that resemble 
(statistically) pure random numbers, with the exception of the 
functional form of the invariant density and the existence of de- 
caying time correlations. The results presented here show that 
the system, in its critical state, feels the details of the underlying 
dynamics, even if preserving the universality class. 

The time correlations in the updating produce a change in the 
non-universal features. In particular we showed that, as these 
correlations increase, the critical state of the system moves to- 
wards a more ordered configuration, that is the threshold is 



higher. This correspondence is made evident, for example, in 
the case of the shifted Tent map. There, a completely deter- 
ministic system reproduces the original BS results. 

We would like to draw the attention of the reader to the com- 
plementarity of the results presented here and those obtained by 
Bianucci et. al. |32| Q . They showed that if a variable w (say a 
Brownian particle) is weakly coupled to a system, provided this 
system is chaotic or ergodic, the resulting deterministic motion 
of the variable w conforms to a standard fluctuation-dissipation 
process. In fact, the irregularities of the deterministic statis- 
tics are washed out by the time scale separation between the 
system of interest (represented by w) and the chaotic subsys- 
tem. The chaotic system is referred to as a "booster" [p2|-p4[. 
This is completely analogous to what happens in the BS model. 
Noise (thermal or otherwise) can be replaced by a determinis- 
tic system without significant changes in the stationary state. 
Stochasticity in the updating rule is sufficient but not neces- 
sary: SOC persists, even in the absence of chaos, for (ergodic) 
periodic updating rules, if in a different universality class. More- 
over, the conditions required from a deterministic system to be 
an appropriate booster are very similar to those required (from 
the updating rule) for SOC to be preserved. 

Summarizing, the results presented here indicate that the fea- 
ture ensuring SOC in systems with extremal dynamics, is not 
the randomness of the updating but the fact that the choice of 
the site where the change is performed (namely the minimum 
rule) is random. Moreover, as long as there is enough diversity 
among the species on the lattice, the longer the memory (or the 
internal correlation) of each member, the higher the threshold. 
Indeed, in the case of chaotic maps, the diversity is ensured by 
the random assignation of the initial values and as much as the 
chaoticity is increased we see that the threshold decreases. In 
the case of the periodic map instead, the random initial con- 
ditions do not provide enough diversity. Indeed, in order to 
have SOC, the internal time-scales, i.e. the periods, have to be 
distributed in a disordered fashion. Briefly, one needs enough 
diversity for SOC to appear jSl]]. 

At this stage, several questions arise. On the one hand, the 
behavior of the threshold with the parameter r in the case of 
Bernoulli updating needs to be explained. Second, and perhaps 
more important from a theoretical point of view, what is the 
exact relation between correlation time and the position of the 
threshold? 

Finally, we believe that these results add strength to the rel- 
evance of SOC in physics and biology, since they allow different 
microscopic mechanisms to underlie its appearance as a collec- 
tive behavior. 
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